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LONG-TERM  GOALS 

The  calibration  and  validation  of  surface  wind  and  stress  retrievals  from  oceanic  synthetic  aperture 
radar  (SAR)  imagery  is  especially  difficult  in  tropical  cyclone  (TC)  conditions.  The  geophysical  model 
functions  (GMFs)  that  characterize  the  radar  backscatter  in  terms  of  the  near-surface  wind  vector  for 
different  viewing  geometries  are  currently  poorly  characterized  for  the  very  high  wind  and  strong 
ocean  surface  wave  conditions  that  are  present  in  all  TCs.  Our  long-term  goal  is  to  develop  a  novel 
methodology  to  use  surface  pressure  observations  and  a  planetary  boundary  (PBL)  model  for 
calibration  and  validation  (Cal/Val)  of  SAR  GMFs  in  TC  conditions  and  to  produce  scene-wide  wind 
vector  retrievals  that  are  most  consistent  with  the  image  backscatter,  the  GMF  and  the  PBL  model. 

OBJECTIVES 

The  objectives  of  this  research  are  to  (1)  develop  the  methodology  for  deriving  TC  SEP  fields  from 
first-guess  surface  wind  vector  estimates  based  on  various  GMF  formulations;  (2)  Use  these  SEP  fields 
to  derive  surface  wind  vector  retrievals  that  are,  in  a  least-squares  sense,  a  scene- wide  optimal  surface 
wind  retrieval  that  is  consistent  with  the  observable  linear  features  in  the  SAR  image,  the  GMF  and  the 
PBL  model;  (3)  Develop  an  optimization  scheme  that  seeks  the  minimum  adjustment  to  the  surface 
wind  vector  field  that  will  minimize  the  difference  between  observed  (e.g.  via  drop  sondes  or  buoys) 
and  SAR-derived  bulk  pressure  gradients  across  the  image.  The  optimized  surface  wind  field  can  then 
be  used  either  to  assess  or  adjust  the  GMF. 
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APPROACH 


GMFs  describe  the  average  radar  backscatter  from  the  sea  surface  in  terms  of  the  surface  wind  speed 
and  the  viewing  geometry  (incidence  and  azimuth  angles).  Scatterometers  obtain  multiple 
approximately  simultaneous  looks  at  multiple  different  viewing  geometries  (at  least  fore  and  aft  looks) 
at  the  same  patch  of  the  sea  surface.  The  incidence  angles  are  known,  so  it  is  possible  to  estimate  the 
most  likely  wind  speed  and  azimuth  angle  (i.e.  wind  direction)  that  best  explains  to  the  separately 
measured  radar  backscatters  from  the  different  viewing  geometries. 

As  compared  to  scatterometer  wind  vector  retrievals,  the  standard  method  for  obtaining  surface  wind 
vectors  from  SAR  imagery  is  very  difficult.  The  same  GMFs  are  used,  but  there  is  only  a  single 
viewing  geometry  for  each  wind  vector  cell.  Wind  directions  are  estimated  by  identifying  linear 
features  in  the  imagery  at  scales  ranging  from  -0.5  km  to  2-4  km  wavelengths.  These  linear  features 
are  assumed  to  be  the  surface  imprint  of  PBL  coherent  structures  that  roughly  align  with  the  surface 
wind  (Foster,  1997;  2005;  Morrison  et  ah,  2005;  Lorsolo  et  al.  2008;  Zhang  et  al.  2008;  Ellis  and 
Businger,  2010).  These  linear  features  cannot  be  identified  throughout  the  entire  SAR  scene  (which  is 
commonly  about  450  by  450  km),  so  wind  directions  are  interpolated  between  the  identified  linear 
features  using  different  methods.  These  directions  and  the  backscatter  are  input  to  the  GMF  and  wind 
speeds  are  retrieved. 

Our  research  is  performed  in  collaboration  with  Chris  Wackerman  of  General  Dynamics  (GD)  and 
Jochen  Horstmann  of  NATO  Undersea  Research  Centre  (NURC).  GD  and  NURC  have  developed 
separate  methods  for  estimating  wind  directions.  The  GD  and  NURC  wind  directions  have  been 
merged  into  a  single  wind  direction  product  and  incorporated  into  a  new  processing  system  (WiSAR) 
that  produces  surface  wind  vector  maps  at  1  km  spacing.  WiSAR  has  been  installed  and  is  being  run  by 
Mike  Caruso  at  the  Center  for  Southeastern  Tropical  Advanced  Remote  Sensing  (CSTARS,  Hans 
Graber  director).  We  use  the  WiSAR  wind  vector  fields  as  input  to  our  SEP  pattern  and  SEP-filtered 
surface  wind  retrieval  system. 

Because  the  pressure  gradient  force  is  a  dominant  term  in  the  PBE  momentum  budget,  the  imprint  of 
the  surface  stress  field  can  be  used  to  estimate  the  surface  pressure  gradient  field  through  the  use  of  a 
diagnostic  PBE  model.  In  its  simplest  form,  PBE  model  assumes  that  the  mean  advective  forces  are 
relatively  small  and  that  the  flow  is  neutrally-stratified  and  baro tropic.  Adding  complexity,  the  PBE 
model  includes  the  effects  of  thermal  winds,  boundary  layer  stratification,  a  gradient  wind  correction 
for  curved  flow  and  momentum  entrainment  across  the  boundary  layer  top.  We  significantly  improved 
the  nonlinear  dynamics  representation  in  terms  of  a  modified  gradient  wind  correction  for  TC 
conditions,  which  added  a  second  tuning  parameter  to  the  PBE  model.  We  also  developed  and  included 
an  important  correction  for  the  inherent  radial  dependence  of  boundary  layer  depth  in  a  strong  swirling 
vortex  consistent  with  the  vortex  boundary  layer  model  that  was  developed  in  Eoster  (2009). 

Once  the  pressure  gradient  fields  have  been  obtained  from  the  surface  winds  and  the  PBE  model,  a 
least-squares  optimization  technique  is  used  to  find  the  best-fit,  zero-mean  surface  pressure  pattern  that 
matches  the  pressure  gradient  field.  If  pressure  observations  are  available,  the  average  difference 
between  them  and  the  zero-mean  field  is  the  least-squares  optimal  estimate  of  what  is  effectively  the 
integration  constant  that  results  from  converting  pressure  gradients  into  a  pressure  field. 

The  premise  of  using  pressure  data  for  SAR  wind  CalW  al  is  based  on  the  following  fact:  Even  without 
using  ancillary  data  to  set  the  absolute  value  of  the  pressure  field,  the  bulk  pressure  gradient  (BPG) 
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between  any  two  points  in  the  SAR-derived  pressure  field  is  the  optimal  estimate  of  that  pressure 
difference  derived  from  the  SAR  imagery.  Hence  pressure  differences  are  more  useful  than  point-by¬ 
point  pressure  comparisons  assessing  the  quality  of  surface  wind  retrievals.  Furthermore,  the  surface 
wind  field  is  driven  in  large  part  by  the  pressure  gradients,  so  wind  vector  errors  and  pressure 
difference  errors  are  directly  related. 

The  derived  SLP  patterns  may  be  used  as  inputs  to  the  PBL  model  to  re-derive  an  “SLP-filtered” 
surface  wind  field.  This  product  is  a  scene-wide  estimate  of  the  surface  wind  vectors  that  enforces 
consistency  between  the  all  of  the  wind  vectors  and  the  pressure  fields.  We  have  extensively 
documented  this  overall  methodology  and  the  high  quality  of  our  SLP  fields  and  derived  wind  vectors 
using  QuikSCAT  and  ASCAT  scatterometer  wind  vector  data  in  several  papers  (Patoux  et  al.  2003; 
2008;  2011). 


Our  applications  to  TC  conditions  with  SAR  wind  vectors  have  shown  that  both  the  directions  and 
speed  in  the  SLP-filtered  winds  can  be  significantly  changed.  In  particular,  the  GMFs  in  high  winds 
work  best  when  the  surface  winds  are  across  the  radar  beam  and  work  poorly  when  the  radar  beam 
looks  either  up-  or  down-wind.  The  SLP-filtered  winds  “fill-in”  the  low  wind  holes  near  the  TC  core  in 
these  up-  and  down-wind  regions.  They  can  also  fill  in  regions  in  the  SAR  WiSAR  product  that  have 
been  masked  due  to  rainy  conditions  or  for  being  out  of  range  of  the  GMFs.  Prior  to  the  ITOP  field 
campaign,  we  had  only  a  handful  of  good  SAR  cases  to  study  with  near-in-time  in  situ  data  such  as 
dropsondes  and  stepped- frequency  microwave  radiometer  (SFMR)  surface  wind  speeds.  Examination 
of  these  cases  has  shown  that  the  SLP-filtered  wind  speeds  are  a  major  improvement  over  the  input 
wind  speeds.  A  qualitative  inspection  of  the  wind  directions  suggests  that  they  are  more  realistic  than 
the  input  directions,  which  tend  to  have  too  little  inflow  in  the  inner  core  region.  However,  validation 
of  wind  directions  is  currently  on-going  using  data  from  the  ITOP  experiment. 


We  now  describe  how  a  SAR-derived  pressure  field  can  be  used  to  calibrate  the  surface  winds.  We 
pose  the  CalWal  problem  as  an  optimization  problem  in  which  we  seek  the  minimum  adjustment  to  the 
surface  vector  field  that  will  minimize  the  difference  between  the  SAR  BPG  and  that  derived  from 
observations.  A  first-guess  surface  wind  vector  field  from  SAR  is  produced  by  WiSAR  from  which  we 
calculate  a  first-guess  surface  pressure  field.  Given  multiple  surface  pressure  observations  within  this 
field,  e.g.  from  dropsondes  or  aircraft-deployed  buoys,  we  can  calculate  multiple,  M  =  A(A-1)/2  , 


BPGs  within  the  image.  The  cost  function  to  be  minimized  is  the  RMS  difference  between  the  SAR- 
derived  BPGs  and  the  observations  along  with  constraints  on  the  corrections  to  the  wind  field  so  that 
the  optimization  procedure  does  not  impose  unreasonable  wind  vector  corrections.  Assuming  M  BPG 
observations  and  K  total  surface  wind  vectors,  the  cost  function  is: 


ivi  i. 
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The  first  term  represents  the  RMS  difference  between  the  SAR  and  observed  BPGs.  The  second  term 
minimizes  the  difference  between  the  SAR  and  SFMR  surface  wind  speeds.  The  third  term  measures 
the  RMS  vector  magnitude  of  the  corrections  to  the  original  surface  wind  field  and  acts  as  a 
smoothness  constraint  -  it  penalizes  large,  localized  wind  vector  modifications. 


Each  optimization  step  adjusts  the  surface  wind  vectors  and  recalculates  the  surface  pressure  field  and 
the  BPGs.  The  optimization  treats  the  U  and  V  components  of  the  SAR  winds  separately  since  their 
error  characteristics  are  approximately  independent  and  Gaussian  (Ereilich,  1997).  The  computation 
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cost  is  nontrivial  because  eaeh  step  in  the  optimization  requires  a  new  nonlinear  surface  pressure  field 
in  order  to  evaluate  the  cost  function. 

It  should  be  emphasized  that  sinee  surfaee  pressure  fields  are  integral  properties  of  the  surface  vector 
wind  fields,  wind  adjustments  must  oecur  over  a  broad  spatial  region  rather  than  just  loeally  near  the 
pressure  data.  Thus,  even  though  we  use  point  wise  pressure  data,  they  imply  wind  correetions  at  a 
large  number  of  wind  veetor  eells. 

WORK  COMPLETED 

Completed  tasks:  (1)  Working  version  of  SLP  retrieval  code,  ineluding  neeessary  PBL  model 
developments,  that  is  eompatible  with  GD,  NURC  and  WiSAR  file  formats  (as  well  as  for  ultra-high 
resolution  QuickSCAT  scatterometer  files)  has  been  developed,  extensively  tested  and  installed  at 
CSTARS.  (2)  SLP  eode  has  been  installed  at  NURC  and  we  have  been  experimenting  with  Horstmann 
to  determine  if  it  ean  (or  should)  be  ineluded  as  an  integrated  part  of  the  NURC  SAR  wind  retrieval 
system.  (3)  Test  version  of  the  CalWal  optimization  scheme  has  been  developed  in  FORTRAN.  At 
present  it  uses  an  older  version  of  the  PBL  model  and  SLP  retrieval  scheme.  Only  preliminary  tests 
have  been  performed.  (4)  Data  from  the  ITOP  field  program  have  been  proeessed  through  the  system 
and  wind  and  SLP  fields  have  been  placed  on  the  ITOP  data  catalog. 

RESULTS 

Much  of  the  researeh  involved  a  colleetion  of  historical  Canadian  RadarSAT-1  SAR  images  from  the 
Hurricane  Wateh  Program  provided  through  an  announcement  of  opportunity  by  the  Canadian  Spaee 
Agency,  NOAA,  NASA  and  CSTARS.  From  this  collection,  we  examined  a  number  of  typhoons  and 
hurricanes.  The  methodologies  were  developed  and  improved  using  these  data  prior  to  the  ITOP  field 
program.  For  the  purposes  of  this  report,  we  foeus  on  the  ITOP  data  from  three  images  of  Western 
pacifie  typhoons  that  span  the  range  from  a  weak  Cat- 1  typhoon  to  a  Cat-5  super-typhoon  and  from 
exeellent  in  situ  data  to  no  in  situ  data. 

Typhoon  Malakas,  22  Sep,  2010,  20:30  UTC 

At  the  time  this  image  was  aequired,  Malakas  had  just  reaehed  typhoon  strength.  It  also  has  the  best  in 
situ  data  coverage  of  any  of  the  images.  The  C-130  aireraft  was  within  the  SAR  scene  at  the  time  of 
aequisition  for  very  eomplete  storm  survey  between  20:00  to  01:00  (23'^‘*)  UTC  that  ineluded  31  drop 
sondes  and  SFMR  winds.  Beeause  the  storm  was  a  weak  Cat-1,  the  wind  speed  range  was  within  reaeh 
of  the  standard  GMF,  CMOD5N.  The  image  was  acquired  in  HH  polarization,  so  we  had  to  apply  the 
Thompson  and  Beal  (2000)  polarization  eorrection  using  the  parameter  a  =  0.8.  The  storm  eenter  was 
in  the  low  incidenee  angle  region  of  the  image,  so  saturation  was  an  issue. 

The  C-130  data  were  mapped  to  the  SAR  image  using  a  two-step  proeess.  First,  we  used  the  HRD 
flight-level  storm  center  estimates  to  find  the  best-estimate  eye  loeation  for  the  eomplete  C-130  storm 
survey.  Similarly  we  ealeulated  the  eye  location  for  each  time  along  the  drop  sondes’  deseents.  From 
these,  we  calculated  the  distanee  and  azimuth  from  the  storm  eenter  for  eaeh  individual  SFMR  or 
sonde  measurement.  An  initial  estimate  of  the  storm  eenter  at  the  SAR  overpass  time  was  interpolated 
from  the  HRD  best  traek.  Adjusted  latitudes  and  longitudes  relative  to  this  eenter  were  ealeulated  from 
the  instantaneous  distances  and  azimuths.  The  HRD  center  estimate  at  the  overpass  time  can  differ 
from  the  center  direetly  estimated  from  the  SAR  image  due  to  errors  in  the  flight  level  estimates,  tilt  in 
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the  storm  with  height  and  the  possibility  of  multiple  eirculation  eenters.  So,  a  final  eonstant  translation 
to  be  applied  to  all  of  the  adjusted  data  loeations  is  obtained  from  the  differenee  between  the  SAR  and 
HRD  eenters.  The  C-130  flight  track  normalized  to  the  SAR  overpass  time  is  shown  in  Figure  1. 

Figure  2a  shows  the  input  wind  field  using  the  merged  GD/NURC  wind  directions  and  the  CMOD5N 
GMF.  This  is  a  descending  pass  and  the  low  incident  angles  are  along  the  right  edge  of  the  image.  The 
winds  were  retrieved  on  1  km  pixels.  The  white  vectors  are  spaced  30  km  apart  and  have  uniform 
length.  The  superposed  black  vectors  are  10  m  wind  directions  from  the  C-130  drop  sondes.  The  sonde 
surface  wind  speeds  are  plotted  as  white-outlined  boxes  with  the  same  colormap  as  the  SAR  wind 
speeds;  they  should  fade  into  the  background  if  the  speeds  agree.  Figure  2c  shows  the  derived  SLP 
pattern  normalized  using  the  drop  sonde  splash  pressures,  which  are  shown  a  color  shaded  white- 
outlined  boxes.  The  SLP-filtered  winds  are  shown  in  Figure  2b.  Note  the  significant  changes  in  both 
wind  speed  and  direction  between  the  SLP-filtered  winds  and  the  raw  CMOD5N  input  winds. 

The  SAR  wind  speeds  are  compared  with  the  SFMR  surface  wind  speed  estimate  in  Figure  3  a.  SFMR 
surface  wind  speeds  are  calibrated  to  1-min  average  winds;  SAR  wind  speeds  represent  10-min 
averages.  The  SAR  winds  are  multiplied  by  1.2  (Harper  et  al.  2010)  to  compensate  for  the  differences 
in  GMF  calibration.  Both  the  raw  and  SLP-filtered  speeds  compare  well,  but  the  SLP-filtered  winds 
track  SFMR  better.  A  first  check  on  the  SAR-derived  SLP  is  shown  in  Figure  3b.  The  C-130  flight- 
level  data  include  an  estimate  of  the  surface  pressure  below  the  aircraft.  The  SAR  SLP  underneath  the 
C-130  tracks  the  C-130  surface  pressures  quite  well. 

The  best  test  of  the  methodology  is  shown  in  Figure  4a.  The  SAR  SLP  pattern  is  normalized  to  the 
drop  sondes  and  the  relatively  small  scatter  shows  that  the  shape  of  the  pressure  pattern  is  good.  A 
better  test  is  shown  in  Figure  4b,  which  compares  the  pressure  differences  between  all  pairs  of  sonde 
surface  pressures  with  the  corresponding  SAR  SLP  pressure  differences.  If  the  shape  of  the  retrieved 
pressure  pattern  is  accurate,  the  slope  of  the  best  fit  line  should  be  unity.  The  in  situ  data  confirm  that 
the  SAR-derived  SLP  pattern  agrees  well  with  the  drop  sonde  data. 

Overall,  the  comparison  between  the  SAR  SLP  field  and  the  SLP-filtered  winds  with  the  C-130  data 
shows  that  the  methodology  has  worked  very  well  in  this  case.  There  are  evident  issues  with  the  wind 
retrievals  in  the  low  incidence  angles  (northeast  side  of  the  storm).  In  these  cases  the  CMOD5N  model 
could  not  match  the  backscatter  and  selected  the  closest  match  to  the  backscatter  for  the  given  viewing 
geometry.  Because  the  SLP-filtering  acts  as  a  scene-wide  retrieval,  it  finds  that  slightly  different  wind 
directions  are  a  better  overall  fit  to  the  inputs.  The  GMF  in  high  winds  and  low  incidence  angles  is 
very  sensitive  to  the  wind  direction.  If  the  SLP  method  is  viable,  we  can  use  the  SLP-derived  wind 
directions  to  improve  the  CMOD5N  wind  retrievals.  We  use  the  SLP-filtered  directions  as  inputs  to 
CMOD5N  and  produce  a  new  wind  field  iteratively.  The  results  are  shown  in  Figure  5.  Figure  5  a 
shows  the  raw  input  winds.  The  dark  blue  patch  to  the  northeast  of  the  storm  center  is  the  region  of  no 
CMOD5N  retrieval.  Figure  5b  shows  the  result  of  a  single  iteration  using  SLP-filtered  directions  and 
Figure  5c  shows  the  result  of  a  second  iteration.  Note  that  the  region  of  no  solution  shrinks  after  each 
iteration. 

Typhoon  Megi  15  October,  2010  21:00  UTC 

Typhoon  Megi  was  a  very  small,  very  strong  typhoon  that  propagated  rapidly  across  the  Pacific.  The 
compact  size  presents  a  major  challenge  to  the  PBL  model  used  in  the  SLP  retrieval,  and  the  results 
show  the  further  research  is  needed  to  improve  its  performance  for  storms  of  this  type.  At  the  time  the 
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image  was  acquired,  the  C-130  aircraft  was  dispatched  from  Guam  to  deploy  a  string  of  buoys  across 
the  predicted  path  of  the  storm.  It  made  a  mini  survey  of  the  storm  on  its  return  flight  about  six  hours 
after  the  satellite  overpass.  So,  there  is  much  less  in  situ  data  and  that  data  have  a  large  difference  with 
the  overpass.  However,  according  to  the  JMA  best  track,  the  storm  strength  did  not  increase 
significantly  during  this  time.  The  aircraft  data  were  adjusted  to  the  SAR  image  time  as  described 
above;  the  adjusted  track  is  shown  in  Figure  6. 

The  input  and  SLP-filtered  surface  winds  are  shown  in  Figure  7ab  and  the  SLP  field  is  shown  in  Figure 
7c.  The  small  size  of  Megi  relative  to  Malakas  is  evident.  The  SFMR  wind  speeds  are  plotted  on 
Figure  7ab  using  the  same  colormap  as  the  SAR  winds,  and  they  evidently  agree  better  with  the  SLP- 
filtered  winds  than  with  the  raw  CMOD5N  winds.  However,  the  sonde  surface  wind  and  splash 
pressures  have  larger  disagreement  with  the  SAR  winds  and  SLP. 

The  SFMR  surface  winds  and  C-130  flight-level  surface  pressure  estimates  are  compared  to  the  raw 
and  SLP-filtered  SAR  winds  in  Figure  8a.  The  agreement  is  best  for  the  SLP-filtered  winds  on  the  first 
SW-to-NE  transect  of  the  eye  and  on  the  final  eye-to-SE  transect.  The  SAR/C-130  surface  pressure 
comparison  in  Figure  8b  suggests  that  the  PBE  model  in  the  pressure  retrieval  system  is  having 
difficulty  capturing  the  very  strong  pressure  gradients  in  this  very  compact  storm.  Some  of  this  error 
might  be  associated  with  a  low  bias  in  the  raw  input  winds.  The  scatter  plots  of  pressure  differences 
between  the  sonde  splash  pressures  and  the  SAR  SLP  are  consistent  with  too  weak  pressure  gradients 
in  this  case  (Figure  9).  Improving  the  model’s  performance  in  compact  storms  is  a  focus  of  the 
research  during  the  follow-on  grant. 

Typhoon  Megi  17  October,  2010  21:41  UTC 

This  image  captured  Megi  near  the  time  that  it  reached  super-typhoon  status  and  just  before  it  made 
landfall  on  Luzon.  JMA  estimated  a  central  pressure  of  885  mb  and  peak  winds  of  64  m/s.  A  C-130 
drop  sonde  from  about  9  hours  earlier  measured  a  peak  near-surface  wind  speed  of  85  m/s.  This  flight 
was  the  closest  in  time  to  the  satellite  overpass,  but  the  JMA  best  track  showed  that  the  storm  was 
continuing  to  deepen  by  at  least  12  mb.  Hence,  for  this  case  we  do  not  have  reliable  in  situ  data  for 
comparison  even  as  it  represents  the  most  severe  possible  challenge  to  wind  and  pressure  retrieval  due 
to  its  small  diameter  and  extreme  strength. 

A  major  issue  with  the  high  wind  behavior  of  C-band  GMFs  in  high  winds  is  they  become  multi¬ 
valued,  with  high  and  low  estimates  of  wind  speed.  Usually  we  use  the  lower  wind  speed  ambiguity.  In 
this  case  however,  it  was  not  possible  to  generate  a  reasonably  deep  SLP  field  from  this  image.  So,  a 
small  number  (fewer  than  200  pixels)  of  high  wind  ambiguities  replaced  the  lower  ambiguities  in  the 
inner  core  region.  The  criteria  were  that  the  radius  from  the  eye  had  to  be  less  than  35  km  and  the  high 
ambiguity  could  not  be  more  than  60  m/s  stronger  than  the  low  ambiguity. 

The  input  winds,  SLP-fdtered  winds  and  SLP  are  shown  in  Figure  10.  The  pressure  field  was 
normalized  using  the  sonde  splash  pressures  from  the  C-130  flight  from  9  hours  earlier  without 
correction  for  further  storm  intensification.  The  most  interesting  feature  of  this  case  is  the  suggestion 
of  a  double  eyewall  in  the  SLP-filtered  winds.  Validating  this  structure  is  the  focus  of  current  research 
in  the  follow-on  grant. 
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Partial  Summary  of  Drop  Sonde  Surface  Pressure  Results 

Figure  1 1  show  a  summary  from  four  Atlantic  hurricanes  and  two  cases  from  ITOP  (excluding  Megi) 
of  the  drop  sonde  splash  pressure  differences  and  SAR-derived  pressure  differences.  Overall  the 
agreement  is  quite  good,  although  further  research  should  reduce  the  scatter.  Because  the  Megi  results 
are  inconsistent  with  these  other  cases  further  research  is  needed  to  improve  the  PBL  model 
performance  for  the  condition  of  small  and  intense  storms. 

Organized  Surface  Stress  Convergence  Lines 

This  research  did  find  strong,  but  preliminary,  evidence  a  new  level  of  organization  in  the  near-surface 
boundary  layer  with  characteristic  wavelengths  of  ~  10  km.  Previous  research  had  documented  that 
low  aspect  ratio  (wavelength/PBL  depth)  PBL-spanning  roll  vortices  are  a  very  common  feature  of 
hurricane  and  typhoon  boundary  layers.  Foster  (2005)  presented  a  theoretical  explanation  why  they 
ought  to  be  the  expected  state  of  the  tropical  cyclone  boundary  layer.  These  rolls  form  at  small  angles 
to  the  mean  PBL  wind  and  have  typical  wavelengths  of  about  2.5  to  4  times  the  PBL  depth,  or 
nominally  1500  to  2500  m.  Hence,  their  effects  should  not  be  resolved  in  these  1  -km  pixel  wind  speed 
images.  (Note  however  that  their  imprint  can  be  resolved  in  the  raw  10  to  25  m  SAR  pixels.  This 
imprint  is  used  to  make  the  first  guess  at  wind  directions  in  the  wind  speed  retrieval  process.)  Although 
the  strongest  direct  effect  of  the  rolls  is  a  periodic  enhancement  and  reduction  of  the  mean  azimuthal 
winds  in  the  radial  direction,  there  is  an  associated  overturning  circulation  that  induces  alternating 
bands  of  convergence  and  divergence,  with  the  same  1.5  to  2.5  km  spacing  as  the  rolls,  in  the  surface 
wind  field. 

One  major  benefit  of  the  scene-wide,  SLP-filtered  wind  retrievals  is  that  it  is  possible  to  make  sensible 
calculations  of  wind  vector  gradients.  Figure  12  shows  surface  wind  convergence  fields  calculated 
from  SAR  images  of  six  hurricanes  and  typhoons.  They  all  show  a  high  degree  of  organization  on 
multi-km  (-'10  km)  scales.  These  patterns  are  roughly  aligned  along  the  surface  wind  direction  and 
hence  bear  a  lot  of  similarity  to  the  PBL  rolls  seen  in  observations  (Wurman  and  Winslow  Morrison  et 
ah,  2005;  Lorsolo  et  al.  2008;  Zhang  et  al.  2008;  Ellis  and  Businger,  2010)  and  theoretically  derived  in 
Foster  (2005).  What  remains  to  be  explained  is  their  very  large  aspect  ratio,  which  is  on  the  order  of  5 
to  10.  Even  though  the  observations  in  Morrison  et  al.  (2005)  were  dominated  by  the  low  aspect  ratio 
roll  modes,  their  observations  documented  that  a  small  fraction  of  their  roll  observations  were  from 
high  aspect  ratio  rolls.  Possible  explanations  are  coupling  with  tropospheric  gravity  waves,  nonlinear 
roll-roll  interactions  and  coupling  to  convection.  Their  dynamical  role  is  also  unknown.  They  may 
contribute  to  nonlocal  (i.e.  non-gradient)  fluxes  of  enthalpy  and  momentum  across  the  PBE  or  they 
might  precondition  the  PBE  and  lower  troposphere  for  future  convection  bands. 

We  recently  became  aware  (Prof.  Gregory  Tripoli,  pers.  comm.)  of  some  numerical  modeling  studies 
that  find  10  km  scale  roll  vortices  throughout  their  boundary  layer  outside  of  the  radius  of  maximum 
wind.  Near-surface  convergence  plots  from  their  numerical  model  are  a  good  match  to  the  SAR 
observations.  They  attribute  these  large  aspect  ratio  rolls  to  a  mechanism  similar  to  that  described  in 
Poster  (2005). 

It  should  be  noted  that  the  theory  in  Poster  (2005)  does  indeed  predict  that  nonlinearly  stable,  large 
aspect  ratio  rolls  (-10)  are  possible.  Pigure  13  reproduces  Pigures  4  and  7  from  Poster  (2005).  The 
interpretation  of  Pigure  13a  is  that  a  wide  range  of  wavenumbers  with  a  range  of  orientation  angles  are 
unstable  (positive  growth  rate).  All  of  these  unstable  modes  are  roll  vortices  that  vary  primarily  in 
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wavelength  and  orientation  relative  to  the  mean  wind.  The  fastest  growing  modes  agree  with  the  low 
aspect  ratio  rolls  that  have  been  documented  in  the  observational  studies.  However,  large  aspect  ratio 
rolls  (corresponding  to  wavenumbers,  a,  near  0.25  are  also  linearly  unstable,  although  more  slowly 
growing  than  the  low  aspect  ratio  roll  modes.  Furthermore,  Figure  13b  shows  that  all  linearly  unstable 
modes  are  nonlinearly  stable  -  that  is,  the  exponential  growth  is  checked  by  self-modification  and 
modifications  of  the  mean  flow  in  a  manner  that  the  rolls  and  mean  flow  equilibrate  into  a  new  mean 
state  with  finite  non-local  fluxes  of  heat  and  momentum.  So,  large  aspect  ratio  rolls  that  could  match 
the  structures  seen  in  the  SAR  in  the  convergence  fields  are  consistent  with  theory  of  Foster  (2005). 

Because  the  large  aspect  ratio  roll  modes  grow  more  slowly  than  the  basic  roll  modes,  all  things  being 
equal,  they  would  not  at  first  be  expected  to  survive  the  “competition”  for  extracting  mean  flow 
energy.  However,  Mourad  and  Brown  (1990)  demonstrated  that  nonlinear  roll-roll  interactions  are  a 
viable  mechanism  for  coupling  energy  from  the  initially  faster  growing  low  aspect  ratio  rolls  into  the 
large  aspect  ratio  rolls.  The  coupled  system  finds  a  new  equilibrium  in  which  the  high  aspect  ratio  rolls 
coexist  with  the  low  aspect  ratio  rolls.  This  mechanism  was  used  to  explain  the  observations  of  cold  air 
outbreaks  off  the  pack  ice  in  which  both  large  and  small  aspect  ratio  rolls  were  found  to  coexist.  A 
similar  mechanism  may  be  in  play  in  the  hurricane  boundary  layer.  This  is  an  exciting  avenue  for  new 
research. 

These  features  are  an  unexpected  result  of  this  research  and  future  research  efforts  will  focus  on 
developing  a  dynamical  explanation  for  their  formation  and  maintenance  and  on  uncovering  what  role 
they  play  in  typhoon  physics. 

IMPACT/APPLICATIONS 

The  strong  correlation  between  the  derived  SLP  and  the  observations  suggests  that  the  CaWal 
optimization  strategy  has  strong  potential.  Those  results  should  have  a  significant  positive  effect  on 
future  SAR  wind  retrievals  in  TC  conditions.  This  will  be  pursued  in  future  research. 

PUBLICATIONS 

We  are  currently  writing  papers  describing  the  methodology  and  results  described  in  this  report. 
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Foster,  RC,  J  Patoux  and  RA  Brown,  “Applications  of  OVW-Derived  Surface  Pressures  in  Tropical 
Cyclones”,  2010  International  Ocean  Vector  Winds  Meeting,  Barcelona,  Spain,  18-20  May,  2010. 

Foster,  RC  :  “An  Innovative  Method  for  Validating  Surface  Winds  Derived  from  SAR  Imagery  for 
Tropical  Cyclones”  RADARSAT  Hurricane  Applications  Workshop,  Canadian  Space  Agency,  John  H. 
Chapman  Space  Centre,  Saint-Hubert,  Quebec,  5  Oct.  2010 

Foster,  RC  and  J  Patoux:  “Using  Surface  Pressure  to  Improve  SAR  High  Wind  Retrieval  in  Tropical 
Cyclones”,  PO25J-04,  Ocean  Sciences  Meeting,  Portland,  OR,  21-25  Feb,  2010. 

Foster,  RC  and  J  Patoux:  “Diagnosing  tropical  cyclone  surface  wind  fields  from  SAR”,  17* 
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FIGURES 


Figure  1:  C-130  storm  survey  of  typhoon  Malakas,  22  Sep, 
2010,  20:00  to  01:00  (12*'^)  UTC  superposed  on  the  SAR 
image  acquired  22  Sep,  2010,  20:30  UTC.  The  track  has 
been  adjusted  to  the  overpass  time  as  described  in  the 
text.  Color  shading  shows  increasing  time  (blue  to  red), 
squares  are  30  m.  apart.  HRD/Willoughby  storm  center 
track  is  shown  in  black;  dots  are  60  min.  apart.  SAR 
storm  center  is  magenta  star;  HRD  center  at  20:30  is 
black  star. 
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Figure  2:  Malakas  22  Sep  2010,  20:30  UTC.  (a)  Raw  input  winds  derived  from  merged  GD/NURC  wind 
directions  and  CMOD5N.  (b)  SLP-filtered  surface  winds,  (c)  SAR-derived  sea-level  pressure  (SLP).  The 
wind  vectors  have  uniform  length  and  show  the  SAR  directions  at  30  km  spacing.  The  black  vectors 
show  the  sonde  10  m  directions.  In  (a)  and  (b)  the  sonde  10  m  wind  speeds  are  shown  in  white  outlined 
squares.  In  (c)  the  boxes  are  the  sonde  splash  pressures. 


11 


Figure  3:  Malakas,  22  Sep  2010  20:30  UTC.  (a)  Comparison  of  SAR  surface  wind  speeds  with  C-130 
SFMR  wind  speeds  (biack).  The  raw  SAR  wind  speeds  are  in  biue  and  the  SLP-fiitered  wind  speeds  are 
in  red.  The  SFMR  rain  rate  is  in  cyan.  The  SAR  winds  have  been  muitipiied  by  1.2  to  correct  for  the 
different  mean  wind  scaiing  (~10  min  for  SAR  and  ~1-min  for  SFMR).  (b)  Estimate  of  surface  pressure 
from  C-130  fiight  ievei  data  (biack).  SAR  SLP  in  red.  Ciosest  sonde  spiash  pressures  to  C-130  sub¬ 
aircraft  point  are  piotted  as  dots  coior  coded  as  distance  in  km  of  the  spiash  point  from  the  aircraft. 

GWCfac  =  0.5200;  RMSE  =  1 .7  mb  GWCfac  =  0.5200 


Figure  4:  Maiakas  22  Sep  2010,  20:30  UTC.  (a)  Comparison  of  sonde  spiash  pressure  with  SAR -derived 
pressure,  (b)  Comparison  of  buik  surface  pressure  differences  between  aii  pairs  of  drop  sondes  and  the 
SAR  SLP. 
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CM0D5:  WiSAR 


CM0D5:SLP  AZ1 


CM0D5:  SLP  AZ2 


Figure  5:  Malakas  22  Sep  2010,  20:30  UTC.  Effect  of  using  SLP-fiitered  wind  directions  to  improve 
CMOD5N  retrievais  in  high  wind  and  iow  incidence  angie  conditions,  (a)  raw  speeds  using  GD/NURC 
directions,  (b)  Retrieved  speeds  using  SLP-fiitered  wind  directions  (1  iteration),  (c)  Retrieved  speeds 
using  SLP-fiitered  wind  directions  (2  iterations).  The  region  of  no  soiution  is  the  dark  biue.  Note  how  it 
shrinks  when  SLP-fiitered  directions  are  used. 


Figure  6:  Typhoon  Megi  15  October  2010, 
21:00  UTC.  Adjusted  C-130  fiight  track. 
Symbois  as  in  Figure  1. 
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Figure  7:  Typhoon  Megi,  15  Oct  2010,  21:00  UTC.  As  in  Figure  2.  Aiso  piotted  is  the  SFMR  wind  speed 
using  the  same  coior  map  as  the  SAR  wind  speed. 
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Figure  8:  Typhoon  Megi,  15  Oct  2010  21:00  UTC.  As  in  Figure  3. 
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Figure  9:  Typhoon  Megi,  15  Oct  2010  21:00  UTC.  As  in  Figure  4. 


BPG  Observed  (mb) 


2010-10-17  21:42:17 
SAR  Input  (m  s’") 


2010-10-17  21:42:17, 
SLP-Filtered  (m  s’^) 


2010-10-17  21:42:17 
Sea-Level  Pressure  (mb) 


Figure  10:  Supertyphoon  Megi,  1 7  October  21 :41  UTC.  As  in  Figure  2. 
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Figure  1 1 :  As  in  Figure  4.  Summary  of  drop  sonde  comparisons  for  four  Atiantic  hurricanes  and  two 
iTOP  typhoons. 
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Multi-km-Scale  Surface  Wind  Conv./Div.  Patterns 
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Figure  12:  Surface  wind  convergence  fieids  caicuiated  from  the  SLP-fiitered  winds  for  six  hurricanes 
and  typhoons.  The  spacing  between  the  bands  is  about  10  km. 


Figure  13:  (a)  Figure  4  and  (b)  Figure  7  from  Foster  (2005).  The  originai  captions  are:  FiG.  4.  Linear 
instabiiity  growth  rate  ag  as  a  function  of  angie  from  azimuthai  wind  e’  =  e-  90°  and  wavenumber  a  for  r  = 
40km  and  AT  =  -  5  K.  Soiid  iines  denote  positive  growth  rates,  dash-dotted  iines  denote  zero  growth, 
and  dashed  iines  denote  negative  growth  rates.  Contour  intervai  is  0.005  (nondimensionai).  The 
equivaient  piot  for  neutrai  stratification  is  very  simiiar  with  overaii  reduced  growth  rates  and  the  peak 
shifted  toward  iarger  e’.  FiG.  7.  First  Landau  coefficient  ai  corresponding  to  the  conditions  in  Fig.  4  for 
ao  >  0.  The  ag  =  0  contour  is  shown  as  dash-dot.  Contour  intervai  is  0.0005  (nondimensionai). 
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14.  ABSTRACT 

Synthetic  aperture  radar  (SAR)  images  of  the  sea  surface  underneath  tropical  cyclones  have  the  potential  to  provide  air-sea 
interaction  data  at  km-scale  resolution.  However,  many  challenges  remain  in  extracting  such  data  in  these  extreme  environments. 
This  research  focused  on  the  problem  of  extracting  surface  vector  wind  and  sea-level  pressure  fields  from  the  SAR  images.  The 
geophysical  model  functions  (GMFs)  that  describe  the  radar  backscatter  are  currently  poorly  characterized  for  the  very  high  winds 
and  sea  states  in  tropical  cyclones.  Calibration  and  validation  of  GMFs  is  severely  inhibited  by  the  lack  of  accurate  in  situ  surface 
wind  data.  We  demonstrate  that  surface  pressure  data  provides  a  more  useful  method  for  validating  and  improving  surface  vector 
wind  retrievals  and  that  scene-wide  wind  vector  retrievals  that  employ  our  technique  are  superior  to  the  standard  pixel-by-pixel 
methodologies  that  are  currently  standard  practice.  Aircraft  and  in  situ  data  from  the  Impact  of  Typhoons  on  the  Pacific  (ITOP) 
2010  field  program  were  used  to  develop  and  demonstrate  the  methodology. 
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